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We present a new formulation of jet quenching in perturbative QCD beyond the eikonal ap- 
proximation. Multiple scattering in the medium is modelled through infra-red-continued (2 — > 2) 
scattering matrix elements in QCD and the parton shower describing further emissions. The inter- 
play between these processes is arranged in terms of a formation time constraint such that coherent 
emissions can be treated consistently. Emerging partons are hadronised by the Lund string model. 
Based on this picture wc obtain a good description of the nuclear modification factor at Rhic and 
Lhc and predict its behaviour at very large px ■ 



High-momentum transfer processes, such as single in- 
clusive hadron spectra, jet-hke particle correlations and 
reconstructed jets, have been advocated for a long time as 
precision tools for studying the dense QCD matter pro- 
duced in ultra-relativistic nucleus-nucleus collisions. The 
rates and distributions of such processes can be bench- 
marked in the absence of medium effects both experimen- 
tally (e.g. through their measurements in proton-proton 
collisions) and theoretically (due to their perturbative 
calculability). In nucleus-nucleus collisions at Rhic and 
at the Lhc, they show generic and strong medium modifi- 
cations commonly referred to as jet quenching. This mo- 
tivates using such hard processes as well-calibrated hard 
probes of dense QCD matter. A central challenge for the- 
oretical work is to relate, in a quantitatively controlled 
and unambiguous way, the observed medium modifica- 
tions to fundamental properties of the QCD matter that 
induces them. In this context, applying a perturbative 
description plays a central role. 

The first measurements of jet quenching signatures at 
the Lhc sparked an intense debate on the extent to which 
a purely perturbative description of jet quenching is fea- 
sible. On the one hand, it seems plausible that several 
features of partonic in-medium propagation may resist 
a perturbative formulation. This is the case, for in- 
stance, for effects of colour flow between jet and medium 
that modify the non-perturbative late-time dynamics of 
jet hadronisation. On the other hand, several perturba- 
tivcily motivated jet quenching models [7 10] resting on 
the eikonal limit [1-6] have been shown to reproduce cen- 
tral features of the nuclear modification factors at Rhic 
and at the Lhc. However, the eikonal limit is known to 
be plagued by large numerical uncertainties [11], and this 
complicates an assessment of whether a purely perturba- 
tive description can account for jet quenching. 

Here, we formulate a conceptually new framework 
based on our understanding of perturbative QCD in non- 
eikonal kinematics and to document its first implemen- 
tation in the medium-modified final state parton shower 
Jewel. This formulation treats the scattering of an en- 
ergetic parton off a medium constituent as a standard 
2 — 2 partonic scattering process described by QCD ma- 



trix elements and valid in the full perturbative regime. 
Further following the standard description of all hard 
scattering events at colliders, large logarithms associated 
to coUinear singularities in real emission matrix elements 
are resummed by the parton shower. This provides a 
systematically improvable approximation to extra gluon 
emissions that also includes radiation off the scattering 
centre and generates naturally configurations with any 
number of additional patrons. The parton shower thus 
introduces radiative energy loss, while the matrix ele- 
ments model the coUisional one. The ambiguity between 
elastic and inelastic scattering is thus resolved, and the 
distinction between vacuum and medium induced radia- 
tion becomes obsolete. Treating all emissions fundamen- 
tally in the same way introduces in a natural way the 
interplay of radiation from different sources. 

In the following we discuss how this idea can be imple- 
mented into a Monte Carlo event generator allowing for 
a democratic treatment of all partonic fragments, and, 
at the same time, accounting in a local and probabilistic 
formulation not only for exact four momentum conser- 
vation including all recoils, but also with known accu- 
racy [12, 13] for the dominant quantum interference effect 
in medium-induced gluon radiation, the so-called non- 
Abelian Landau-Pomerantchuk-Migdal (Lpm) effect. 

The parton shower evolution is governed by the Su- 

dakov form factor, which can be interpreted as the prob- 
ability for having no emission between two scales Q-^ and 
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The functions Pba{z) are the (lunegularised) Altarelli- 
Parisi splitting functions describing the energy fraction 
z taken by parton h in the splitting a b -\- c. The 
coupling is evaluated at the transverse momentum of the 
splitting; it is given by fc^. — ~ for initial state 
and by ~ z(l — z)q^ for final state emissions. 

In order to include secondary scatters and initial state 
radiation off such collisions, we introduce 'partonic pdf 's' 
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fj(x, Q^), which analogously to the proton pdf's describe 
- at leading order - the probability that a parton j with 
energy fraction x has been radiated from a parton i when 
the latter is resolved at a scale Q^. These pdf's are 
also employed in the typical backward evolution of initial 
state radiation, as known, e.g., from simulating proton- 
proton collisions. Since the ccntre-of-mass energies of 
scatterings in the medium are typically rather small, it is 
unlikely that more than one initial state splitting is initi- 
ated in secondary scatters. In this approximation [14] the 
partonic pdf's can be integrated analytically, for instance 
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J^qq(a:) , (2) 



and analogously for the other partonic pdf's. 

The interaction of the parton shower with the parton 
constituents in the medium depends on properties of the 
medium. We specify the cross sections for 2 — >■ 2 pro- 
cesses as 



a,iE,T)= J d\i\ / da: 



je{q,q,g} 



dt 



(3) 

where the pdf takes into account possible initial state 
radiation off the energetic projectile. Note that here, 
implicitly, we neglected a similar evolution experienced 
by the target, i.e. there is no initial state radiation off 
the medium parton. The maximum momentum transfer 
|i|max is determined by the kinematics of the scattering. 
Neglecting the scattering centre's momentum one finds 
l^lmax — 2ms(T)[i?p — nip], where ms{T) stands for the 
(temperature dependent) scattering centre's mass and Ep 
and Trip are the projectile parton's energy and mass, re- 
spectively. The boundaries on the x-integral are obtained 
from the requirement that k'j_ > Qq/4: and are given by 
Xmin{\i\) = Qo/(4|t|) and x^^^{\i\) = 1 - Qo/{Mi\)- For 
the partonic cross section we keep leading terms only, but 
we regularise them with a Debye mass /xd ~ 3T. They 
therefore read 

da-,,,-,, ^ TT 9,,-, 9 , + (s — 1^1)^ 

-^(..m = c.^.im+A) ^^^^^ 

Each parton emitted in subsequent parton showering 
has a finite formation time that is parametrically given by 
r « w/A;^ i=a {x)E/Q'^. When the formation time is larger 
than the mean free path in the medium, emissions from 
all scattering centres within in the formation time inter- 
fere destructively giving rise to the non-Abelian Landau- 
Pomcrantchuk-Migdal (Lpm) effect. This has been in- 
corporated following the algorithm of [12, 13], with the 



only difference that in the current framework gluon for- 
mation times cannot overlap. Compared to the original 
algorithm, this amounts to a mild and well understood 
suppression of induced radiation compared to the ana- 
lytic result without changing the shapes of the gluon dis- 
tributions [15]. 

When emissions from different sources, i.e. initiated 
by different hard scatterings (where one may be the 'vac- 
uum' parton shower) compete with each other, it is solely 
the gluon with the shortest formation time that will be 
emitted. This criterion is enough to regulate the inter- 
play of all (vacuum and medium induced) emissions. 

To obtain the example results below, the following set- 
up has been used: In all cases, the final state parton 
shower is generated by Jewel [16], where elastic scat- 
tering and medium-induced radiation is implemented as 
discussed above. Initial dijet production is simiilated us- 
ing the matrix elements and initial state parton shower 
of Pythia 6.4 [17], the latter running in a virtuality- 
ordered mode to provide the best fit to the Jewel sim- 
ulation. In all results, the strong coupling is running at 
one loop with as(mz) = 0.128, consistent with findings of 
other leading order calculations. The infra-red cut-off of 
the parton shower was chosen as Qo = 1-5 GeV to yield a 
good description of Lep data. The Cteq6l1 [18] pdf's 
as provided by Lhapdf [19] are used; simulating nuclear 
collisions the nuclear modification of Eps09lo [20] is em- 
ployed. The medium is modelled using a simple variant of 
the Bjorken model [21], which roughly accounts for the 
main features of the produced matter. It describes the 
boost-invariant longitudinal expansion of an ideal quark- 
gluon-gas. The transverse profile is chosen such that the 
energy density is proportional to the density of wounded 
nucleons. Jewel can in principle be interfaced with any 
medium model. For the hadronisation the strings are 
built by Jewel and then handed over to the Pythia 
hadronisation routine [22, 23]. The colour flow and con- 
sequently the strings can be constructed in different ways. 
The default choice reported here is to treat the struck 
scattering centres as if they were gluon emissions. The 
scattering centre is then always colour-connected to the 
projectile, before the parton shower sets in. 

In principle, the modelling of interactions in the 
medium is free of parameters, as the parton shower and 
its infra-red cut-off Qq, hadronisation and the couplings 
are fixed from e+e~ and pp collisions. In practice, how- 
ever, the infra-red regulator in the cross section (4) adds 
an additional freedom. Here, we choose /Ud ~ 3T, in 
accordance with parametric estimates. The uncertain- 
ties related to this choice will be discussed in detail else- 
where, but they do not affect the conclusions drawn in 
the present study. 

In the absence of a medium the Jewel parton shower 

reduces to a standard (vacuum) parton shower. The 
parton shower implementation in Jewel was thoroughly 
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tested in e+e~ collisions at Lep and pp reactions at Rhic 
and Lhc. As an example for the generally good per- 
formance, Fig. 1 shows a comparison with the neutral 
pion spectrum measured by Phenix [24], which forms the 
baseline for the measurement of the nuclear modification 
factor i?AA- Above p_L ~ 4GeV, the Jewel+Pythia 
results agree with the data on a level of roughly 15 % 
over about 6 decades. 
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FIG. 1. Neutral pion spectrum in pp collisions at ^/s = 
200 GeV simulated with Jewel+Pythia and compared to 
Phenix data [24]. The Analysis of Monte Carlo events and 
plots were made with Rivet [25]. 



Turning to medium-modifications of hadron spectra, 
we fix the critical temperature Tc — 165 MeV, consistent 
with the expected temperature at which QCD matter 
transfers from partonic to hadronic degrees of freedom. 
To obtain a fair agreement with the measured nuclear 
modification factor at Rhic, Jewel-|-Pythia requires 
an initial temperature of Ti = 350 MeV at initial proper 
time Tj — 0.8 fm, see Fig. 2. This is remarkably con- 
sistent with the input parameters in fluid dynamic sim- 
ulations of heavy ion collisions. We note that at high 
p±, where the Monte Carlo results are reliable, they re- 
produce both the factor ^ 5 suppression and the ap- 
proximately flat pj^-dependence seen in the data. Vary- 
ing the Debye mass by ±10% has a significant effect on 
the overall suppression, indicated by the grey band, but 
hardly affects the shape. Other sources of uncertainties 
like pdf-uncertainties, choice of initial time Tj etc. can be 
expected to be smaller and will be discussed elsewhere. 

The charged hadron multiplicities measured in heavy 
ion collisions constrain the initial entropy density of the 
system, SiTi oc dN/dy, Si oc ei/Tj oc 7]'^ and therefore 
allow to relate the initial temperatures at Rhic and at 



PHENIX data, 0-10% centrality 
JEWEL+PYTHIA 



0.6 - 



0.4 



0.2 



■ ■ ■ - ■ ■ ■ 















1 




i 





10 15 
P, (GeVl 



20 



25 



FIG. 2. Nuclear modification factor for neutral pions in 
Au-I-Au collisions at ^ = 200 A GeV in the 0-10% central- 
ity class simulated with Jewel-|-Pythia and compared to 
Phenix data [26] (only statistical errors shown). The grey 
band indicates a variation of the Debye mass by ±10%. 



the Lhc, 
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The observation of a factor 2.2 increase in the event mul- 
tiplicity from Rhic to Lhc is therefore consistent with 
an initial temperature Ti = 530 MeV at ri = 0.5 fm 
at the Lhc. There is some freedom in initializing the 
fluid dynamic evolution at the Lhc at a different initial 
time Ti, but this is numerically unimportant. At early 
times the parton shower is dominated by emissions at 
rather high scales initiated by the initial hard scatter- 
ing, and this high virtuality protects the partons from 
medium-induced emissions and makes them insensitive 
to the medium at early times. Thus, the medium at the 
Lhc is specified in terms of parameters fixed in Fig. 2. As 
seen from Fig. 3, the calculation of Jewel±Pythia then 
leads without any additional adjustments to a very good 
agreement with preliminary data of the nuclear modifi- 
cation factor at the Lhc. 

To understand the characteristically different p±- 
dependencies of Raa at Rhic and at the Lhc, we have 
performed control simulations in which a Lnc-like dis- 
tribution of hard processes is fragmented in a Rnic-like 
medium, and vice versa (data not shown). This showed 
that the significant change in the slope of Raa{p±) from 
Rhic to the Lhc can be attributed fully to the y/s- 
dependence of the distribution of initial hard processes. 
We therefore conclude that a purely perturbative dy- 
namics of parton energy loss supplemented by an ar- 
guably simple model of the medium whose characterisa- 
tion matches physical expectations, can account for the 
main features of the measured nuclear modification fac- 
tors, including the strength of the suppression pattern. 
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CMS preliminary data, 0-5% centrality 
ALICE preliminary data, 0-5% centrality 
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FIG. 3. Nuclear modification factor for charged hadrons in 
Pb+Pb collisions at ^ = 2.76 A TeV in the 0-5% centrality 
class simulated with Jewel -|-Pythia and compared to pre- 
liminary Cms [27] and Alice [28] data (only maximum of 
statistical and systematic and errors shown). 



and its ^/s- and pj_-dependence. 

At very large pj_ the nuclear modification factor contin- 
ues to rise above unity. This is a purely kinematical effect 
that becomes visible at very large p±^ where the energy 
loss starts to vanish. The elastic scattering of energetic 
partons converts longitudinal momentum into transverse 
momentum and thus effectively makes the pj^-spectrum 
harder. While this is a generic effect its size and the 
turn-on point to some degree dependend on the medium 
model, as they are sensitive to the amount of scattering 
centres encountered in the forward direction. 

In this publication we have presented a novel descrip- 
tion of jet quenching dynamics entirely based on stan- 
dard perturbative technology also used in the simulation 
of proton-proton collisions, and implemented it into the 
Monte Carlo generator Jewel. In contrast to other sim- 
ilar attempts, we do not only overcome certain limita- 
tions imposed by the eikonal limit [29], but we also re- 
solve the dichotomy in the description of in-medium and 
vacuum radiation, basing all emissions on the same par- 
ton shower. We reinterpret induced radiation as regular, 
although infra-red-regulated, 2 — >■ 2 parton scatterings, 
supplemented with our parton shower. This of course 
makes any distinction of elastic and inelastic scattering 
obsolete. We hinted at various ways of how our model 
can systematically be improved, using a well-understood 
perturbative language. This framework is flexible enough 
to accommodate observables such as correlations, or jets, 
which reach far beyond the very inclusive ones studied 
here. For the fairly successful simulation of the nuclear 
modification factor at Rhic and Lhc we used a simple 
Bjorken model. Taking the description of this inclusive 
data as a proof of principle, we reserve further and more 
detailed investigations to future work. 
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